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ABSTRACT 

When a cold shell sweeps up the ambient medium, a forward shock and a reverse shock 
will form. We analyze the reverse-forward shocks in a wind environment, including their 
dynamics and emission. An early afterglow is emitted from the shocked shell, e.g., an optical 
flash may emerge. The reverse shock behaves differently in two approximations: relativistic 
and Newtonian cases, which depend on the parameters, e.g., the initial Lorentz factor of the 
ejecta. If the initial Lorentz factor is much less than IIAeI^'^ Aq ^* the early reverse 
shock is Newtonian. This may take place for the wider of a two-component jet, an orphan 
afterglow caused by a low initial Lorentz factor, and so on. The synchrotron self absorption 
effect is significant especially for the Newtonian reverse shock case, since the absorption 
frequency Va is larger than the cooling frequency v,. and the minimum synchrotron frequency 
Vra for typical parameters. For the optical to X-ray band, the flux is nearly unchanged with 
time during the early period, which may be a diagnostic for the low initial Lorentz factor of 
the ejecta in a wind environment. We also investigate the early light curves with different wind 
densities, and compare them with these in the ISM model. 

Key words: shock waves - gamma rays: bursts - stars: wind 



1 INTRODUCTION 

Long-duration gamma-ray bursts (GRBs) may originate from the collapse of massive stars (Woosley 1993; Paczynski 1998). The probable 
associations between GRBs and supernovae have been detected in several cases, e.g., the most confirmed ORB 980425 / SN 1998bw (Galama 
et al. 1998, Kulkarni et al. 1998) and ORB 030329 / SN 2003dh (Hjorth et al. 2003), which give a firm link to the coUapsar model. The light 
curves of several afterglows also show the supernova component, such as ORB 970228 (Reichart 1999), 980326 (Bloom et al. 1999), 011121 
(Bloom et al. 2002, Greiner et al. 2003), 021004 (Schaefer et al. 2003), and so on. Thus, the surrounding environment is wind-type. Much 
work about wind-type environment analyses has been done (Dai & Lu 1998; Meszaros, Rees & Wijers 1998; Panaitescu & Kumar 2000; 
Chevalier & Li 2000), including the features of the afterglow light curves and the comparison with the interstellar medium (ISM) model. 

At the beginning of the interaction between the ejected shell and environment, an early afterglow will emerge from the reverse shock, as 
predicted by Meszaros & Rees (1997) and Sari & Piran (1999b). The prompt optical emission from GRB 990123 (Akerlof et al. 1999) and 
GRB 021211 (Fox et al. 2003; Li et al. 2003) were observed, though observations of a very early afterglow were difficult before the Swift's 
launch (Gehrels et al. 2004). The optical flash of GRB 990123 was immediately analyzed (Meszaros & Rees 1999; Sari & Piran 1999a). 
They pointed out that the optical flashes mainly come from the contribution of the reverse shock. Then, plenty of theoretical analyses were 
advanced. The dynamics, numerical results of optical and radio emission, and analytical light curves, from the reverse and forward shocks in 
uniform environments were discussed (Sari & Piran 1995; Kobayashi & Sari 2000; Kobayashi 2000). Early afterglows in wind environments 
were also considered by several groups (Chevalier & Li 2000; Wu et al. 2003; Kobayashi & Zhang 2003). Panaitescu & Kumar (2004) 
considered the reverse-forward shock scenario and the wind bubble scenario for the two observed optical flashes from GRB 990123 and 
GRB 021211. However, little discussion comes to a Newtonian reverse shock in wind environment. We consider this case in the following 
sections. 

Kobayashi, Meszaros & Zhang (2004) noticed the synchrotron self-absorption (SSA) effect on the early afterglow. We find that other 
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parameters like the thickness of the shell and the initial isotropic kinetic energy can also influence the self absorption, even up to the optical 
wavelength. In this paper, we derive the complete scaling-laws of the SSA frequency for all cases. 

In the simulations of Zhang, Woosley & MacFadyen (2003), the initial Lorentz factor can be as low as about tens. For the structured 
jet model (Kumar & Granot 2003), it is likely that the jet has low Lorentz factors at the wings of the jet. Huang et al. (2002) considered 
that a jet with an initial Lorentz factor less than 50 may cause an orphan afterglow. Rhoads (2003) also pointed out that the fireball with a 
low initial Lorentz factor will produce a detectable afterglow, though no gamma-ray emission is detectable. There are indications that some 
GRBs ejecta have two components: a narrow ultra-relativistic inner core, and a wide mildly relativistic outer wing (Berger et al. 2003; Huang 
et al. 2004; Wu et al. 2005; Peng et al. 2005). When the mildly relativistic shell collides with the wind environment, the reverse shock is 
Newtonian. Kobayashi (2000) considered the Newtonian reverse shock in a uniform environment. However, no systematic analysis has come 
into the Newtonian reverse shock in a wind environment. In this work, we discuss the Newtonian reverse shock, which is mainly caused by a 
low initial Lorentz factor. In this case, the optical emission flux from the Newtonian shocked region exceeds that from the relativistic forward 
shocked region. 

Some authors have used the early afterglow as a diagnostic tool of gamma ray bursts' parameters for the ISM case (Zhang, Kobayashi 
& Meszaros 2003; Nakar et al. 2004) and for the wind case (Fan et al. 2004; Fan et al. 2005). Accordingly, the behavior of early reverse- 
forward shocks should be completely described for the wind case. We derive the analytical scaling-laws of dynamics and radiations for both 
relativistic and Newtonian reverse shock cases in the wind environment in §^|2|and|3] and give the numerical results of radio to X-ray light 
curves in ^ We present some discussions in ^ 



2 HYDRODYNAMICS 

Let's consider a uniform and cold relativistic coasting shell with isotropic kinetic energy Eo, Lorentz factor 74 = + 1 ^ 1, and width in 
observer's frame Ao, ejected from the progenitor of the ORB. This shell sweeps up a free wind environment with number density ni — Ar^^, 
where 77 is the initial ratio of Eo to the rest mass of the ejecta (Piran, Shemi & Narayan 1993). The interaction between the shell and the 
wind develops a forward shock propagating into the wind and a reverse shock propagating into the shell. The two shocks separate the system 
into four regions: (1) the unshocked approximately stationary wind (called region 1 hereafter), (2) the shocked wind (region 2), (3) the 
shocked shell material (region 3), and (4) the unshocked shell material (region 4). By using the shock jump conditions (Blandford & McKee 
1976, BM hereafter) and assuming the equality of pressures and velocities beside the surface of the contact discontinuity, the values of the 
Lorentz factor 7, the pressure p, and the number density n in the shocked regions can be estimated as functions of ni, 714, and r;, where 
714 ~ Ef,/{r)^-Kr'^^4,/^QmpC?) is the comoving number density of region 4. 

Analytical results can be obtained in both relativistic and Newtonian reverse shock limit. These two cases are divided by comparison 
between / and 7I, where / = 714/^1 is the ratio of the number densities between the unshocked shell and the unshocked wind (Sari & Piran 
1995). As shown by Wu et al. (2003) for the wind environment case, / — //(jy^Ao), where I = E()/{A-nA'mp(?) is the Sedov length. If 
/ 74, the reverse shock is Newtonian (NRS), and if / <^ 7I, the reverse shock is relativistic (RRS). 

As discussed by Kobayashi & Sari (2000), even for NRS, the adiabatic index of the post-shocked fluid can be taken as a constant 
7 = 4/3, because the electrons are still relativistic. Then the shock jump conditions can read (BM; Sari & Piran 1995) 

e^ln-irapt? = 72-1, 712/71.1 = 472 + 3, (1) 
e^/n-impC = 73 - 1,713/714 = 473 + 3, (2) 

where rrip is the proton mass, 62 and 63 are the comoving energy densities of region 2 and region 3 respectively, and 712 and 713 are the 
corresponding comoving number densities of particles, which are assumed to consist of protons and electrons. The relative Lorentz factor 
between region 3 and region 4 is 



Assmning 72 = 73, and 72,74 ^ 1, 73 can be expressed as 73 ~ (74/72 + 72/74)/2. The asymptotic solution is 73 ~ -i=7|'''^/^ 

73 ~ ^74/ V"^''* for RRS, while 73 ~ 74 and 73 - 1 ~ f 7!/"^ for NRS. 

The time it takes the reverse shock to cross the shell in the burster's frame is given by (Sari & Piran 1995) 



There are two simple limits involved in the problem: NRS and RRS, in which we can get analytical results. The relative Lorentz factor 73 is 
constant in the whole reverse-shock period for RRS. tA can be derived as tA ~ a/S.o'^Af^^'^ / c, and the corresponding radius of the shell at 
time tA is ta — ctA ~ aAo^Af^^"^ — oiy/TKo, where the coefficient a = 1/2 for RRS and a — i/y/lA for NRS. We will discuss both 
cases separately in the following. 




(3) 




(4) 
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2.1 Relativistic Reverse Shock Case 

In the RRS case, / ^ 74 (i.e., 77 2> 114i553* A^jj^yl, using the relation between the observer's time and the radius ~ (1 + 



z)r/2'y^c, where z is the redshift of the GRB, we obtain T ~ (1 + z)Ao/2c ~ 16.7(1 + 2;)Ao,i2S as the RRS-crossing time in the 
observer's frame. We adopt the conventional denotation Q = Qk x 10*° in this paper except for some special explanations. Using e2 — 63, 

72 =73, together with the above equations, we get the scaling-laws of the hydrodynamic variables for time < T, 

73 ~ 1.9772.5^53 Ao^ia^,;.! (5) 

73 = 72 ^ Sl.5EU'Aoy,''A:]!t (6) 

63 = e2 ^ 1.4 X lO^iS-/^ Ao,?/^a!/!, (I) " erg cm^ (7) 

N,,3 ~ 2.1 X 10'''E53V2.1y (8) 

iV,2 ^ 3.5 X W'^eU'K%AI(1,'-^ (9) 



where A, — 3 x 10"^^ cm~^, and A'^e.i is the number of electrons in the shocked region We note that 73 and 73 do not depend on time. 
This is the property of wind environments, since the densities of the shell and the ambient environment have the same power-law relation 
with radius r (n oc r~^). 

After the reverse shock crosses the shell (t^ > T), the shocked shell can be roughly described by the BM solution (Wu et al. 2003; 
Kobayashi & Zhang 2003; Kobayashi et al. 2004), 

,-3/8 J.-9/8 J.-3/2 ,1/4 -r ,0 /,n\ 

J'i OC ' ,713 (Xt^ ' , 63 OC ' ,r (Xt^ , Ne,3 OC t^, (10) 

,-1/4 ,-5/4 ,-3/2 ,1/2 ,r ,1/2 ,,,, 

72 OC tgg ' , n2 OC ' , e2 OC ' ,r oct^ , Ne,2 oc t^^ . (11) 
These variables can be scaled to the initial values (t® — T), which are given by the expressions for the time t(g < T. 



2.2 Newtonian Reverse Shock Case 

In the NRS case, / 2> 74, the time for the reverse shock crossing the shell is T' ~ tr, — 2.9 x 10^(1 + z)E53r]^'t,A^\-^^s in the observer's 
frame, if we consider the spreading of the cold shell (Piran et al. 1993). The evolution of the hydrodynamic variables before the time T' are 

7-3 - 1 ^ 0.57^ (12) 
73 = 72 ~ 74 (13) 

63 = 62 ~ 5.8E^3^r,l,Al_, (^^j ergcm-3 (14) 

iVe,3 = 2.1 X 10'''i?537?r5 (^^j (15) 

iVe,2 ^ 6.6 X 10^2 £537)1-5^ (16) 

What should be noted is that the values for NRS are not suitable for mildly relativistic reverse shock case. Nakar & Piran (2004) showed the 
difference between the approximated analytical solution and the numerical results in the case of uniform environments. And for the spreading 
of the shell, / decreases with radius. At the crossing time, 73 ~ 1.57 (see equation il2i ). which deviates from the Newtonian reverse shock 
approximation. More accurate values should be calculated numerically. 

After the NRS crosses the shell, the Lorentz factor of the shocked shell can be assumed to be a general power-law relation 73 oc 
(Meszaros & Rees 1999; Kabayashi & Sari 2000). However, the forward shock is still relativistic, and can be described by the BM solution. 
The dynamic behavior is the same as the one in the RRS case. The scaling-law of the two regions are 

73 oct-«/(i+^^',n3 0cC^+«)/^'i+^«', 



63 oc C<«+«'''^(^+^^',r oc t'J^'+''\N.,3 oc t%, (17) 

,-1/4 ,-5/4 

72 OC ' , n2 OC ' , 

62 oc tZ^^^ ,r (X tl{^ , Ne,2 oc . (18) 
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3 EMISSION 

We now consider the synchrotron emission from the shocked material of region 2 and region 3. The shocks accelerate the electrons into 
a power-law distribution: N{"f^)d'y^ = N-^'y~^d'y^{-^^ > 7m), where 7^ is the minimum Lorentz factor of the accelerated electrons. 
Assuming that constant fractions and e_B of the internal energy go into the electrons and the magnetic field, we have B = ^Snesei, where 
Bi is the internal energy density of the shocked material. Regarding that the comoving internal energy of the electrons can also be written as 
eeC = A'^^7~^7emeC^d7e, and the comoving number density n = A'"-y7~^d7e, one can get 7™ = ee(7— l)(mp/me)(p — 2)/(p — 
1), where 7 = 73 or 72 corresponds to the reverse or forward shock, and A'^^ — n{p — l)7m~^- 

The cooling Lorentz factor 7c is defined when the electrons with 7c approximately radiate all their kinetic energy in the dynamical time, 
i.e., (7c — l)meC^ = P{yc)tco, where P{ye) = (4/3)crTc(7e — l)(-B^/87r) (Rybicki & Lightman 1979) is the synchrotron radiation power 
of an electron with Lorentz factor 7e in the magnetic field B, tco is the dynamical time in comoving frame (Sari et al. 1998; Panaitescu & 
Kumar 2000). Then the cooling Lorentz factor 7c — Qivmec/ {(jTB^tco) — 1. 

The synchrotron radiation is taken to be monochromatic, and the corresponding frequency of an electron with Lorentz factor 7e is 
Ve ~ (3/2)7ei^L, where ul ~ qeB / {2-nmec) is the Larmor frequency, and Qe is the electron charge. The critical frequencies are — 
3(1 + 2:)~^77m?eB/(27rmec) and Uc — 3(1 + z)~^'yy'^qeB/{2iTmeC), in the observer's frame respectively, where 7 is the bulk Lorentz 
factor of the emitted region. Before the reverse shock crosses the shell, !^c,2 ~ !^c,3 is satisfied for the two regions having the same energy 
density e, Lorentz factor 7, and the same comoving time tco- 

The synchrotron self-absorption effect should not be ignored, especially at low frequencies, where the emission is modified enormously 
for the large optical depth. Wu et al. (2003) have given the SSA coefficient and the corresponding spectral indices. We here quote the results 
of Wu et al. (2003) and derive the SSA frequency for all six cases in the following. 

The initial distribution of shock-accelerated electrons is 

^(7e) = N^'yl^ 7"! < 7e < 7 

max ■ (19) 

Taking into account the synchrotron radiation energy losses, the power-law distribution of electrons is divided into two segments (Sari et al. 
1998), i.e., 

N{'y,) = N,{ (,+1) (20) 

7e 7e fm, 



for the fast-cooling case (7c < 7m), and 



N{je) = { _(p+l) (21) 

7e ' 7e > 7c, 



for the slow-cooling case (7c > 7m). 

The self-absorption coefficients in different frequency ranges are 



ts 



ci7i~ 








(p- 




C27r 








C372~ 



"1 



-5/3 
-(p+4)/2 
-5/2 



e-"^"^ v->V2 



^^^'^ =1 = 9x2V3r(i/3) p+t/3 - C2 = ^-^2^l\p + f )r(2E±2)r(3£±iO), C3 = 2^^(p + 2), and V2 are the typical synchrotron 
frequencies of electrons of Lorentz factor 71 and 72, and r(a::) is the Gamma function (Wu et al. 2003). 

Electrons in both segments contribute to the SSA. For simplicity, the less important segment is neglected. In general, the absorption 
coefficient is dominated by the electrons between 7c and 7™, for the frequency less than max(i/c, Vm)\ and dominated by the electrons with 
Lorentz factor greater than max(7c,7m), for the frequency larger than max(!/c, Vm)- So, the third expression in equation \21\ is always 
unimportant and can be neglected. We can obtain analytical expressions for the SSA frequency Va by taking kvL — to. 



C2tiV,7r^''^+*Vr+^'^^i ^i«/.„«^2 (23) 

2 

(Pl+4),,(P2+4)/2 L 1 P2+4 



[c2tiV^72 



l^a > 1^2, 



where L = Ne/ (Aivr^n) is the comoving width of the emission region, tq can be defined equal to 0.35 (Frail et al. 2000), 71 — min(7c, 7m), 
72 = max(7c, 7m), Pi is the power-law index of the electron distribution between 71 and 72 (pi = p for slowing cooling, pi — 2 for fast 
cooling), and p2 = p + 1 is the index of the electron distribution with Lorentz factor greater than 72. 

Because the peak spectral power Pi/,max ~ (1 + z)'yme(? ar B / {2iqe) in the observer's frame is independent of 7e, the peak observed 
flux density can be given by F^^max ~ NePu^ma.^/ (^t^D^) at the frequency min(i^c, Vm), where D is the luminosity distance of the gamma- 
ray burst. 
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3.1 Relativistic Reverse Shock Case 

Using the above expressions, we obtain the typical frequencies and the peak flux density in the shocked shell and the shocked wind for the 
RRS case, 



„3 ~ 5.9 X 10'''(l + 2) 



-1-2 1/2 



-1/2 2 



'=£'=3. -1-^53 »?2.5-4,,-lAo^i2 



1/2 j iffi 

r 



Hz, 



!/™,2 1.0 X 10 (1 + z) e,,ej^ _-^E^!^ Aoy^ I 



Hz 



1 c w inl2/i I \-l -3/2 T^l/2 . 1/2 .-2 Iffi 

/^c,3 ==; 1.5 X 10 (1 + z) eg '_^E,Ji ^0,12^,,-! Y ^^2' 
95.3(1 + z)4^_^E,yir^^lAl%A-^^^^D^i Jy, 



F„,2 ^ l.Q{l + z)^il^Eli^A^.,-^Al\'^-'D^^'}y 



-1/2 , 



(24) 

(25) 

(26) 
(27) 
(28) 



where = ee,-o.5 ■ 3(p — 2)/(p — 1). Note that Fi,,max,3 > -F'i^,max.2, i.e. region 3 dominates the emission for the early afterglow, mainly 
because the number of electrons in region 3 is much larger than that in region2. 
We give the scaling-law of the SSA frequency in region 3, 



1^0,3 



1.1 X W^il+z 

2.2 X 10"(1 -Fz 
4.7 X 10"(1 -hz 
2.2 X 10"(1 -Fz 

1.1 X 10^^(1 -fz 

5.2 X 10"(1 -Fz 



-1^6/5 



S,-l-'^53 



l/10„-3/5 .19/10 A -19/10 (t^ . „ 
^2.5 ^.,-1^0,12 \~ ' 



-lpl/6„-l/3 .1/6 A -5/6 / t(B 
-^^53 ^.S ^.,-1^0,12 I ~ 



-2/3 



Hz 



-11/15 



-1^2/5 1/10 pl/30 2/15 .1/3 .-23/30 /48 
^3,-1^53 '/2.5 -1^0, 12 I T 

-111/5 p2/5 -8/5 ,2/5 A "7/5 / £e^\ " " „ 
'^3.-1^53^2.5 -1^0, 12 It) 



-1-6/13^9/26 ^-1/26^2/13 ^9/13 -25/26 ( 



B, -1-^53 '/2.5 -1^0, 12 \~ 

z/o i/lO pl/30 2/15 .1/3 

^3,-1-^^53 '/2.5 -1^0, 12 



Uc < Ua < Vm 
Hz I^c < Vm < fa, 
Va < I^m < Vc 
<Va < Vc 



(29) 



Hz 

-1^-2/5^1/10 ^1/30^2/15^1/3 ^-23/30 ( t^\'^^/^'- V^<V^< Va, 



Here and in the following expressions for Va, we take p — 2.5. If more than one expressions above satisfy the followed restriction, the largest 
i/a is the true value. 



In region 2, 



Va,2 



9.6 X 10"(1 + zr\'^:'^,E^r'Al]l^i^o%' (^) " Hz 



5.7 X W^'\l + zy^Al('^^A~ 



-2/3 / t 



-2/3 



Hz 



-1--1 1/5 



''3,-1-^53 



E. 



-2/5 .6/5 



.D/o A -3/5 / *e 
^*,-l'^0,12 I ~ 



4.5 X 10^^(1 + z) 

1.9 X 10^5(1 + ^)-l,-6/13^|/2_6^^1/26^8/13^-2^^^^ 



7.8 X W'\l + z)-^el^'e]i^!!XrAt('_\A-%' (4^ 



-11/15 
-1 

Hz 

-1 

-11/15 



Hz Vc < Vm < Va, 
Va < V,n < Uc 
Hz I'm < l^a < 
Hz Vm < Vc < Va- 



(30) 



After the reverse shock crosses the shell (t® > T), the behavior of both shocked regions can be described by the BM self-similar 
solution. The power-law indices of emission variables with time are given in table Q Another frequency Vcut should be introduced here 
(Kobayashi 2000) to substitute i^c for no fresh electrons. Vc^t has the same time profile as Vm- If Vm > Vcut, Vm comes down to !/cut for the 
synchrotron cooling. And if Va > Vcut, Va comes down to Vc^t too for no electrons distributed greater than corresponding 7cut. These are 
all represented in columns labeled (S) (for NRS case) and ® in TableQ 

The scaling-law indices of flux densities with time are sophisticated, as they vary with time when any two of v, Ua, Vc{oi i^cut), Vm 
cross each other, where u is the observed frequency. These indices are given in Table |2| For the case v > v^ut, the flux density decreases 
exponentially with observed frequency v, then we take it to be zero, which is denoted by a short horizontal line in Table|2| The numerical 
results will be given in Sj4] 

For the typical parameters, the order of the frequencies at = T is Um > v > fa > i^c for both region 3 and region 2, if the considered 
frequency — 4.55 x lO^^Hz. The flux density from the shocked shell and shocked environment are 



5.4(1 + zf'\l]L\Etn^lA::llA-,%'D^i Jy. 



(31) 



F.., ^ 0.1(1 + zf/\-^L\E%'A-^^,^D-,i Jy. 



(32) 
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3.2 Newtonian Reverse Shock Case 

For NRS case, before the reverse shock crosses the shell (t^ < T'), 



Vm,2 



.max, 2 



4.1 X lQ^\l + zy'eyi'_,E^inl,K[U ^ Hz 



1^4 ^3/2 t®_ 

T' 

1.3x10 (1 + z) eets _i-E'53 77i,5yl,;_i I — 



Hz 



o T inl3/i I \-l -3/2 -2 i-5/2 ^0 it 

Vc,3 = i^c,2 2.7 X 10 (1 + 2) <:s,_i£^53?7i. 5^,-1 ^ Hz 



I ^ 1/2 3 i3/2 r,-2 / i® 

7.6(1 + 2)es;_ir?i.5A.,_i-D28 1^ 
0.2(1 + z)e'J-'_,r^l,Al(l,D2s' Jy 



-1/2 



Jy 



(33) 

(34) 

(35) 

(36) 
(37) 



The SSA frequency in region 3, 
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4^0,3 < 



1 --1^1/5 



-1„4 49/5 7 tfg 
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5.8 X 10«(1 + z)-'e'J''4ZEs^vLAt%f ( \ 



23/iO 
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Hz 
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Hz 

-7/15 



7.9 X 10"(1 + z)-'el/''eX:^^E-,''/'%lTAl]^\'' (^) ' Hz 



Vc <Va < Vm 
Vc < Vm < Va, 
Va <Vm < Vc 
<l'a < Vc 
Vm < Vc < Va, 



and in region 2, 



2.8 X 10"(1 + zr^T-iE^'vlTK^l 
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2.3 X 10^2(1 + z)-ii?-//V5'^.,-i(^ 
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-11/15 
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-1 12/5 ^9/5 f t, 
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[ 1.6xlO-(l + z)-^ei/^^;_°,i5-3"/i%J°f<^\''(^j ^ Hz vm<vc<v., 

After the reverse shock crosses the shell (t^ > T'), the temporal indices of the typical frequencies and the observed flux density are 
also given in tablesQand|2| 

The frequency relations at time = T are v > Uc > Va > Vm for region 3 and Vm > v > Uc > Va for region 2,ifu = 4.55x10"Hz. 
The corresponding optical flux density from region 3 and 2 are 



Hz 

-11/15 



Va < Vc < Vm 
Vc <Va < Vm 
Hz Vc < Vm < Va, 
Va < Vm < Vc 
Vm < Va < Vc 



(38) 



(39) 



TT' cc/i I \-l/4-3/2 1/8 i-,-1/4 5 ^11/8 r-i-2 t 

F^,3 ~ 65(1 -f z) ' ej efl -i-Egg' m.5^,,li-D28 mJy, 



(40) 



n />o/i , n1/2 -1/4 ril/2 il/4 r-i-2 t 

F^,2 ~ 63(1 + z) ' £b^LiE^3 '7i-5^,,_i£>28 mJy- 



(41) 



4 NUMERICAL RESULTS 

The above analytical results can give approximate behaviors of variables as functions of time or frequency, but they are valid only in 
relativistic or Newtonian limits. In the mildly relativistic case, the analytical values deviate from the actual ones very much (Nakar 2004). 
For precise results, a numerical method should be engaged. 

Combining equations Q, J2j and J3j, and using the assumption of equalities between the Lorentz factors and pressures beside the 
surface of the contact discontinuity, one can obtain solutions of 72, 73, 73, 62, 712, 113 numerically. Before the reverse shock crosses the shell. 
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Table 1. The temporal indices for the evolution of v^, v^, Va and -Fi/,max- The notations denote respectively: early, NRS, region 2; (2) early, NRS, region 
3; (D early, RRS, region 2; @ early, RRS, region 3; (5) lately, NRS, region 2; ® lately, NRS, region 3; ® lately, RRS, region 2; ® lately, RRS, region 3. Vc 
is actually Vcut at the columns (6) & ®. 



variable 


t < tA 










notation 




©(D® 


(D 


®® 


® 


® 




Va 


-1 
1 

-2 
2 
3 


1 
1 

23 
10 
5 
6 


3 
2 

1 
2 

8 
5 
2 
3 


15 g+24 

14 9+7 

15 9+24 

14 9+7 
33 9+36 
70 9+35 

15 9+24 

14 9+7 

15 9+24 

14 9+7 
33 9+36 
70 9+35 

(15 9+24) p+32 9+40 
(14 9+7) p+56 9+28 

15 9+24 

14 9+7 
11 9+12 

14 9+7 


15 
8 
15 
8 

_ 3 
5 
15 
8 


Va <.Va<, Vm 
Vc <,Va Vm 


(fast cooling) 

Va 

(slow cooling) 


P+3 
P+5 
-1 

-1 

P+3 
P+5 




p— 6 
P+5 
23 
10 
p-7 
p+4 
p— 6 
P+5 
1 
2 


3 p+5 

2 p+10 

3 
5 

3 p+6 

2 p+8 

3 p+5 
2 p+10 

1 
2 


15 
8 
3 
5 

15 p+26 
8 p+32 

15 

8 

9 

8 


Va Vm <i Va 
Va <.Vm <. Vc 
Vm <^Va^Vc 
Vm <^ Vc Va 



Table 2. The temporal indices for the evolution of flux density (Fi, oc i^). The notations are the same as in tableQ The short horizontal lines indicate the 
radiation vanishes at those cases. 
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3 p''-2 
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1-2 
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3 p-2 
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5 
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1-2 


P^i 








3 p-2 


2 


2 


4 


2 


3 


1 


1 


5 





3 


6 


p-1 


p-2 


3 p-1 


2 


2 


4 


p-2 


p-1 


3 p-2 


2 


2 


4 


2 


3 


1 


5 


5 


7 


2 


2 


4 


p-1 


p-2 


3 p-1 


2 


2 


4 


p-2 


p-1 


3 p-2 


2 


2 


4 


2 


3 


1 


5 


5 


7 


2 


2 


4 


p-2 


p-1 


3 p-2 


2 


2 


4 



< I'a < I'c < i^m 
Va < V < Vc < Vm 
Va < Vc < V < Vm 
Va < Vc < Vm < V 

V < Vc < Va < Vm 
Vc < V < Va < Vm 
Vc < Va < V < Vm 
Vc < Va < Vm < V 

V < Vc < Vm < Va 

Vc < V < Va 
Vc < Vm < Va < V 

V < Va < Vm < Vc 
Va < V < Vm < Vc 
Va < Vm < V < Vc 
Va < Vm < Vc < V 

V < Vm < Va < Vc 
Vm < V < Va < Vc 
Vm < Va < V < Vc 
Vm < Va < Vc < V 

V < Vm < Vc < Va 

Vm < V < Va 
Vm < Vc < Va < V 



5 9+8 
14 9+7 
6 9+4 
14 9+7 



19 9+36 
14 9+7 



5 9+8 
14 9+7 
6 9+4 
14 9+7 
(15 9+24) p+7 9 
28 9+14 

5 9+8 



14 9+7 
25 9+40 
28 9+14 
(15 9+24) p+7 9 
28 9+14 



19 9+36 

14 9+7 
53 9+96 
28 9+14 



21 
8 



1 

2 
_ 1_ 
2 

15 p+3 
16 



1. 
2 

23 
16 
15 p+3 
16 



21 
8 

57 
16 



the value of 73 should be solved from the following equation without approximation, 



(73- 1) (473 + 3) = 



7374 



1 - 



73 V 72 



47374 1 - 



+ 3 



(42) 



and then the other variables can be derived directly. 

We take the parameters 77 = 300, Eo = 1.0 x 10^^ erg, A, = 0.1, Ao = 5.0 x 10^^ cm, ec = 0.3, en = 0.1, and D = 1.0 x 10^** 
cm, for the RRS case. For the NRS case, we set i] = 30, while keeping the same other parameters as in the RRS case. Following the above 
analysis, we can get the emission from the two shocked regions, of which the optical magnitude at frequency v — 4.55 x 10^"* Hz is shown 
in Figure |2| [later] for RRS case and NRS case respectively. The reverse shock dominates the emission at the beginning and fades after the 
shock crosses the shell, which is identical for both RRS and NRS. This effect may be the cause of the so-called optical flash. 
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Figure 1. Flux density at v = 8.46 GHz as function of time. Parameters are rj = 300, Eq = 1.0 X 10^^ erg, A, = 0.1, Ao = 5.0 X 10^^ cm, €e = 
0.3, eg = 0.1, p = 2.5 for the upper panel (RRS case). Only r; = 30 is different for the lower panel (NRS case). The long dashed and short dashed lines 
represent the emission from region 3 and region 2, respectively, and the solid line is the total flux density from both regions. 

After the reverse shock crosses the shell, we choose the parameter (/ = 1 for the dynamic evolution of NRS. Kobayashi and Sari (2000) 
discussed that g should satisfy 3/2<(/<7/2in the ISM environment. A similar conclusion can be drawn for the dynamics of the ejecta in 
the wind environment. As the NRS caimot decrease the velocity of the ejected shell effectively, the shocked ejecta should be quicker than the 
one in RRS case, which satisfies 73 oc r~^^^. On the other hand, the ejecta must lag behind the forward shock, which satisfies 72 oc r^^''^. 
So the range of g should then obey 1/2 < g < 3/2. What's more, the evolutions of the hydrodynamics and the emission do not depend on 
the value of g sensitively. The evolution of 7 with the observer's time has a narrow range from t^^'^ to t'^^'' corresponding to the range of 
9- 

Figures fTl3l show the light curves at radio (8.46 GHz), optical band (4.55 x 10^'* Hz) and X-ray (1.0 x 10^** Hz) respectively. The upper 
panel denotes the RRS case, and the lower panel denotes the NRS case. At low frequencies, Ua is always greater than the observed frequency, 
so the emission at these frequencies is affected by the synchrotron self-absorption enormously, and can be estimated as thermal emission at 
this band (Chevalier & Li 2000). The radio flux density increases with time before and shortly after the crossing time, as shown in FigureQ 
which comes from the increasing number of the accelerated electrons. The flux will be intense enough to be detected if the distance is not so 
large, as the flux is inversely proportional to the square of the luminosity distance. 

The numerical results are well consistent with the analytical ones. For the typical parameters and v = 4.55 x 10^* Hz as the observed 
frequencies, at the crossing time, the orders of the typical frequencies are Vc,i < i^a,3 < 1^ < Vm,,z for RRS case, Vc,'j, < i'm,z < Va,-i < v 
for NRS case, and Vc.i < Va,2 < v < Um,2 for both cases. From Table|2| we find that the corresponding temporal indices are 1/2, —1/4 and 
1/2 for the time before the reverse shock crosses the ejected shell, where p = 2.5. In Figure|2| the slopes can be seen from the four dashed 
lines before the break point, which is the crossing time. The value of the flux density from region 3 at time t = T' is however not consistent 
with the value (55 mjy) given by equation <40> . which is about 3.5 mjy in the figure, since the reverse shock is mildly relativistic. The curves 
are not accordant well with the approximated analytical slopes either. ^ 

For the optical band, the reverse shock dominates the emission at the beginning, and decays quickly after the crossing time, since there 
are no fresh shocked electrons to produce the emission. This is the same for both RRS and NRS cases, as seen in Figure |2| However, the 
X-ray afterglow is always dominated by the forward shock, especially for the NRS case, since the reverse shock is not strong enough, and 

^ The reverse-forward shock is assumed to begin at the fireball's coasting period, which is the initial time for the early-afterglow. However, the coasting radius 
is not zero, though it can be neglected at late times, which has been adopted in the scaling-law analyses. Therefore, at early times, the curve in Figure l2lfor 
RRS case is not straight. As the curve for NRS in the figure begins at 10 s, the influence of the nonzero initial radius can be neglected now. 
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Figure 4. Magnitude as function of time at 4.55 X 10^'* Hz. Parameters are tlie same as in FigurelTlexcept for rj. Tlie T] is 25, 35, and 50 from top to down. 



can't accelerate the electrons to a high stochastic Lorentz factor to emit numerous X-ray band photons. Figure|3|shows the emission at X-ray 
band for both RRS and NRS cases. From these three figures, we can see that the main emission is approximately at optical band. 

As seen in equation <40> . the flux density depends on r) very sensitively. We plot the magnitude as function of time for different -q values 
in FigureO Taking into account the lower panel in Figure|2|and the first two in this figure, we can find that, with other parameters unchanged, 
the larger 77, the larger the flux density, as the flux density is proportional to about rf'^ if p = 2.5. When -q — 50, shown in the lowest panel 
in Figure|4| the reverse shock becomes mildly relativistic. In the relativistic reverse shock case, it is approximately inversely proportional to 
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Figure 5. Flux density in region 2 (solid line) and in region 3 (dashed line) as function of the observed frequency, at the crossing time. Parameters are the same 
as in Figure □with r) = 30. The aiTows on the left of the vertical line denote the time behavior of the con'esponding typical frequencies before the crossing 
time, and the arrows on the right denote the time behavior after the crossing time. 
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Figure 6. Light curves in optical band (4.55 X lO^^Hz) for different wind parameter A. Parameters are the same as in FigureQexcept A, and rj = 30. The 
crosses ( X ) indicate the crossing time of the reverse shock. 

rj. The flux density descends with the increase of rj. Another phenomenon is that, with the increase of the rj, the time for the emission from 
region 2 to overtake the one from region 3 postpones, and then region 3 almost dominate the emission during the whole early period. Because 
the number of the electrons in region 3 is much larger than that in region 2, the emission is dominated by region 3, when the reverse shock 
is powerful enough to accelerate the electrons to emit enough optical band synchrotron photons. Thirdly, no distinct ascending of the optical 
light curves appears before the crossing time for the NRS case. On the contrary, the light curves will descend at the beginning if the emission 
is dominated by the region 3. 

We plot the spectrum at the crossing time for reverse shock and forward shock respectively in Figure |5| The spectrum is a typical 
synchrotron spectrum for the electron energy distribution with index p — 2.5. The breaks are smoothened by the time equal arrival effect 
(Sari 1998). Both curves have three typical frequencies Va, Vc and Vm. The time behavior of the frequencies before and after the crossing time 
is illustrated in this figure. Before the crossing time, Vm,i oc f ""'^ , which is different from the analytical result for NRS case Vm,3 oct^ (listed 
in Table □. This comes from the fact that the term 73 — 1 can't be taken to be much less than 1, especially for a shallow Newtonian reverse 
shock. Consequently, fa/i becomes approximately oc t^'^'^, not oc in TableQ For these parameters, Um,3 is occasionally equal 

to i/c,3 at the crossing time. After the crossing time, as no fresh electrons supply, electrons with stochastic Lorentz factor greater than 7c 
disappears. The maximum electron Lorentz factor jcut varies with time like 7m (Kobayashi 2000). The corresponding Vcut does so. As one 
can see, the maximum typical frequency in region 3 is i^a, so the synchrotron self-absorption effect is important for the reverse shock in NRS 
case. 

The wind parameter A is important for the reverse-forward shock. How do the light curves vary if the wind density varies? We give a set 
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Figure 7. Magnitude at 4.55 X 10^'* Hz as function of time in ISM environment, ni = lcm~^ is the number density of ISM material, rj = 300 for the upper 
panel is RRS case, and = 30 for the lower panel is NRS case. Other parameters are the same as in FigurefTI 

of light curves of the early afterglow with different A in Figure|6| The parameter A is taken from 3 x 10^^ cm~^ to 3 x 10^^ cm~^ {A^ = 1). 
With the increase of A, the reverse shock converts from Newtonian to relativistic. For the extreme NRS case, the emission is dominated by 
the forward shock (the lower light curves in Figure|5l. This makes the light curves (summation of region 2 and region 3) have no break at the 
crossing time. But with the increase of A, the early emission are gradually dominated by the reverse shock, so the breaks (in the upper light 
curves) appear at the crossing time. Another phenomenon is that, as A decreases, the crossing time becomes longer, which is mainly due to 
the spreading of the ejected shell. With the approximation by ignoring the spreading effect, the analytical crossing time for RRS is a constant 
(1 + z)Ao/2c. In Figure|6| we can see that the crossing time converges to (1 + 2)Ao/2c = 167((1 + z)/2)Ao,5xioi2S with the increase 
of A. 

The ISM environment case has been investigated enormously (Sari & Piran 1995; Meszaros & Rees 1997; Kobayashi Sl Sari 2000; 
Kobayashi 2000; Zhang et al. 2003). We here calculate the light curves of the reverse-forward shock for the ISM density n\ = lcm~^ 
and with the same other parameters as the typical values in the wind environment (see Figure 0. For these parameters, the flux densities 
both for the RRS case (upper panel) and for the NRS case (lower panel) increase during the early period, which are different from the wind 
environment case. We can see in Figure|2|that the light curves almost keep unchanged for early times at the optical band. This may be caused 
by the decrease of the number density of the wind. Therefore, whether there exist a rapid increase of the early optical afterglow may be used 
to distinguish between the ISM and wind environments. 



5 CONCLUSIONS AND DISCUSSION 

We have investigated the whole evolution of the dynamics and emission of the reverse-forward shock in a wind environment, by considering 
both the RRS and NRS cases. The temporal indices of the physical quantities are given in TableQand|2| which cover all interrelations of the 
typical frequencies Va^i^m, i^c and the observed frequency v. The flux densities of the emission at radio, optical and X-ray bands as functions 
of time are shown in FiguresQ-|3|by numerical calculations. 

For the ISM model, there exists a transition radius Rn (satisfies /{Rn) ~ 74) for thick shell (Sari & Piran 1995), where the reverse 
shock becomes relativistic. An enormous difference between the wind model and the ISM model is that the Lorentz factor of the shocked 
regions before the crossing time is constant with time for the wind model in RRS case. Since the transition should satisfy 7/74 = 1, the 
ratio / and 74 are both constant, for the wind environment, and thus no transition exists. Therefore, the relative Lorentz factor of the reverse 
shock 73 is also independent of time. Taking into account these properties, we find that the temporal indices are relatively reliable, even 
if the estimates of the Lorentz factors 73)72> and 73 deviate from the actual values, which are caused by the Newtonian and relativistic 
approximations. 

An optical flash emitted from the shocked shell appears for the typical parameters, but perhaps no darkening can be observed at early 
times. A rapid decay occurs after the reverse shock crosses the shell, and then the emission is dominated by the shocked environment 
material. A radio flare lasts for a longer time. It increases continuously even shortly after the crossing time. The X-ray band emission is 
always dominated by region 2. At the optical band and X-ray band, the flux is nearly unchanged at early times especially for the NRS case. 
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which may be used to diagnose the NRS in a wind environment. For the reverse shock, the synchrotron self-absorption can not be neglected, 
since it may exceed the other two frequencies Uc and Vm as the number density increases. 

There is also a possibility that no prompt optical emission is detected. It may be caused by a low initial energy, a low environmental 
density, or strong absorption, and so on. These will decrease the flux density of the early afterglow to go beyond the detector's limits. In the 
5w!/f's era, many early optical and X-ray afterglows will be detected by UVOT and XRT, like GRB050525A (Shao & Dai 2005, Klotz et al. 
2005), and then the parameters may be determined by early afterglow data more precisely. 
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helpful discussions. This work was supported by the National Natural Science Foundation of China (grants 10233010 and I022I001), and 
the Ministry of Science and Technology of China (NKBRSF G 19990754). 
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